NASA-CR-202410 


y 3 ^' ^ -3 


Application of QFT to the Problem of Failed In-Flight 
Controllers During Approach and Landing of a B-720 Aircraft 

Hwei-Lan Chou 
Daniel J. Biezad 

Cal Poly State University 
San Luis Obispo, CA 93407 


Abstract 

Previous studies by NASA Dryden have shown 
the use of throttles for emergency flight control to be 
extremely difficult, especially for landing. Flight control 
using only the throttles to achieve safe landing for a large 
jet transport airplane, the Boeing 720, is investigated using 
Quantitative Feedback Theory. Results are compared to an 
augmented control developed in a previous study. The 
controller corrected unsatisfactory open-loop 
characteristics by increasing system bandwidth and 
damping, but improving the control bandwidth 
substantially proved very difficult. The pitch controller is 
robust in conditions of no or moderate turbulence. The 
roll controller performed well in conditions of no 
turbulence, but is sensitive to moderate turbulence. 
Handling qualities of the augmented control for approach 
and landing were evaluated by piloted simulation flights. 

Notation 

Throttle-Only Flight Control 
Quantitative Feedback Theory 
Velocity-pitch rate derivative 
Sideslip-roll coupling derivative 
pitch rate (deg/sec) 
flight path angle (deg) 
angle of sideslip (deg) 
bank angle (deg) 
pitch angle (deg) 
thrust (lbs) 

natural frequency 
damping ratio 

stick input(full defiection = ] unit) 
pitch rate feedback gain 
flight path angle feedback gain 
sideslip angle feedback gain 
bank angle feedback gain 
transfer functions 
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(a) short form of (s+a) 

1£,cj) short form for s ^ + 2 Qcos + uP 

c.g. center of gravity 


Introduction 

Through throttle manipulations, engine thrust has 
been found useful in providing some controllability for 
multiengine aircraft in emergency situations with severe or 
complete flight control system failures. This paper focuses 
on a particular application, a simulation augmented control 
developed by NASA for Boeing-720 aircraft Throttles-Only 
Flight Control (TOFC) ^ NASA has found the use of 
throttles feasiblcrfor emergency flight control for a range of 
airplanes, and their analyses for a variety of aircraft are 
available in the literature This controller was 
implemented on a high fidelity B720 flight simulator and 
obtained generally good pilot ratings by increasing the bare 
airframe Dutch-roll and phugoid damping. The primary 
aim of this study is to present an alternative control design 
technique based on Quantitative Feedback Theory(QFT) to 
further improve the Dutch-roll damping and to increase 
controller bandwidth for better handling qualities 

The QFT technique ^ was chosen because of the 
insights it provides throughout its design process. It allows 
designers to specify a desired closed-loop frequency 
response with a desired control bandwidth and damping 
characteristic. It shows why the desired performance may 
not be achieved within the given control actuation and rate 
limits. 


In this paper a QFT augmented throttle-only flight 
path controller for approach and landing is presented. 
Complete details of the aircraft model and justification for 
TOFC are not included, but the reader is reminded that 
"trimming* must be possible and "controllability* must 
exist. Augmented control design using QFT is presented in a 
summary fashion. The full justification and step by step 
procedure may be found in Reference 8. 
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B-720 Linear Model 

The empirical transfer function developed for the 
engines is given in short form notation by 

Q x(lbs) 275 

~ (.55X5) 

Severe bandwidth attenuation occurs beyond frequencies of 
1 rad /sec. For this application his prevented the increase the 
closed-loop bandwidth beyond 1 rad /sec within the range of 
available thrust (see Ref. 6). 

Four configuration variations for the B-720 were 
considered as described in the Appendix. They are 
characterized in both the longitudinal and lateral axes by 
excessive resonance / low phase and gain margins, low 
crossover frequency, and large phase angle roll-off. 

QFT Controller Design 

To apply QFT, the aircraft model is rearranged in a 
unit feedback form as shown in Figures 1 and 2. The inner 
pitch rate and sideslip loops were closed using Kq=60 and 
Kp = 4, which were the settings for the original simulation 
augmentation scheme. 


Performance Specification 

QFT allows designers to specify a desired closed - 
lu .^p frequency response with an upper bound B^ a lower 

bound and a tolerance <5 R specified to obtain robust 
performance. The maximum M m is also given to obtain a 
desired system damping. 


Table 1. QFT Performance Specification 


Freq.(r/s) 

0.1 

0.3 

05 

0.7 

1.0 

2.0 

5.0 

Bu(dB) 

17.0 

17.0 

17.3 

-16.0 

-4.0 

20 

-13.0 

BUdB) 

16.8 

-15.0 

12.3 

4.6 

-12.4 

7.1 

-23.0 

5 r (dB) 

0.2 

2 

5 

8.4 

8.4 

9.1 

15.0 


The performance specification shown in Table 1 are 
the desired closed-loop responses for both the y - and <p * 
loops. These two feedback loops are piloted open-loop 
systems. Additional specifications are ususallly given for 
piloted systems, such as a desired control bandwidth of 2 
rad /sec. (see Ref. 9 for transport aircraft landing 
requirements) and a k/s slope near the crossover frequency. 
These added requirements promote good pilot handling 
qualities. 10 
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Design Constraints 


Four configurations were used to study the 
approach and landing of B-720 throttle-only flight control as 
summarized in Table 2 and in the Appendix. Configuration 
1 was used as the nominal confiuration for control design. 


The nominal configurationfcunfig. 1) is. 

c r<*g> confj „ !- .Ul(.203)[.37, 3.01 1 

e.(deg) 8- (.562) [.624, .111] (.441, 157] (5.25) 


Table Z Flight Configurations for B-720 


(Gear Up) 


Conf 

Weight 

(lbs) 

Alt 

(Ft) 

(MSL) 

A/S 

(Kts) 

Flaps 

(%) 

1 

140,000 

4,000 

160 

0 

2 

140,000 

4,000 

145 

30 

3 

160,000 

4,000 

’ 175 

0 

4 

140,000 

4,000 

155 

30 


The performance bounds constraint is a curve on 
the Nichols Chart that shows the performance tolerance, 5 R , 
from Table 1 at each specified frequency. Satisfying the 
tolerance constraint guarantees that the variation of the 
system response due to plant uncertainties will be no greater 

than <5 r There is a performance bound for each frequency. 

U contours are also shown on the Nichols Chart. 
The U contour is a M-circle that has the magnitude of M m , 
with part of the circle stretched for uncertainty at high 
frequencies.-By having the open-loop frequency response not 
penetrate the U contour, the system's damping will be 
guaranteed no less than the damping selected for M m as a 
design constraint. 


and the min. and max. are : 

r r « ko -0053 (.162) [35 , 3.01] 

*-<*•> mm (.40) (.42 , 1.48] (.66 , .01] (5.19) 


maX ' 


.01 (.28) [.46, 3.43] 


(58) (.45, 157] (.92, .14] (5.24) 

Figure 3. Longitudinal Mode Parameter Variation 


Th of the nominal configuration(config. 1) is: 

„♦<<)«) . . .09 [.47 , 3.65] 

C A.<* **> n ° mina = (.98) (.81 . .15] (.26 , 1.07] (5.02) 

and the min. and max. of are : 

/-♦(dc*) • = .06 [.45 , 3.65] 

' (. 98 ) (.60 , .15] [.24 , .93] (5.01) 

c .(« 09 1 61 , 4.33] 

(1.03>|].0, 20][29, 1.091(5.02) 

Figure 4 Lateral Mode Parameter Variation 


Controller Design Technique 


For inner-loop transfer functions G l and G* the 

' "in Pin ' 

parameter variation given by the four configurations can be 
expressed minimums and maximums. There are tradeoffs 
between plant parameter variations and performance. The 
wider the parameter variation, the more restricted the 
constraints; consequently, more compensation is required. 

In this application, due to the engine response limits, the 
performance specification will have to be relaxed because 
there is not enough control power to provide all the 
compensation that is required to meet specifications. 

The minimum and maximum values of transfer 

functions, and G^ , form the uncertainty template. 

The QFT control package ^ allows the designer to input 
maximum and minimum plant parameter variations, but 
due to the software's limitation of handling the quantity of 
uncertainties, some of them were averaged. These variations 
are listed in Figures 3 and 4. 


Poles/zeros/gain compensation will be required to 
reshape the open- loop transfer functions from GJ and 

Go . On the Nichols Chart, gain raises the transfer function 
Pin 

curve, a zero bends the curve to the right, and a pole bends 
the curve to the left. The com pen sat ion -selected forms the 

controller, G c . G^ and G^ after reshaping, become 

respectively and the open-loop transfer functions 
of the flight path angle and bank angle feedback loops. 

Who,. L;, = G e - * and f„ = Gf; * GJ_ . L', 

and Lg p should be kept on and above the Bo(jo>;), for each 
frequency, Wj , on Lg y and to assure robust 

performance. and must also not penetrate the U 

contour in order to obtain the desired damping In this 
application the additional constraint existed which required 
the controller to be physically realizable (zeros not 
outnumbering poles). 
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Lateral Bank Angle Controller 


Longitudinal Flight Path Angle Controller 

Transfer function and its performance bounds, 

Bo(ja>j), and U contour are displayed on a Nichols Chart in 

Figure 5. Since all the frequency points for G^ are below 

their corresponding Bo(ja>j).reshaping is required. Pure gain 
compensation first raises the curve until it touches the U 

contour as shown in Figure 6(G 7 ‘° =16). Compensation is 
then added to avoid the U contour while satisfying all the 
Bo(ja>j) constraints. Note that no realizable poles/zero 
compensation could be found to do this. A zero at .1 rad /sec, 

for example, pulls the whole G^ curve to the right of the U 
contour. But then a pole at any location would make the 
GJ curve penetrate the U contour. Since a compensator 
with only one zero is physically ununrealizable for this 
application, the compensator, G c , for the flight path angle 

feedback loop, G* 1 " , is a pure gain of 16. The performance 
bound, hence the system robustness, was left unsatisfied. 

The frequency response of the close-loop transfer 
function, T v r , where 

» In 

T r ; = L r er / (1 + L r e r ) = K; * G*. ) / (1 + G J * G^ ) 
is shown in Figure 7. It can be seen in Figure 7 that the 
spread between T max and T m j n , has exceeded the <5 R over 

the frequency range .1 to .7 rad /sec. This is because!^ did 

not satisfy the performance bounds over that frequency 
range. The frequency plot of the close-loop response after 
adding a prefilter is shown in Figure 8. A pure gain prefilter 
of 6 proved most effective in increasing the bandwidth and 
in meeting the prescribed specification. 


Transfer function G^ and its performance bounds, 
Botjo^), and U contour were similarly analyzed on a Nichols 
Chart. G* was not only below all the Bo(jcoj) but also 

no 

penetrated the U contour. A controller, 

Gg£ =(s+.l5)/<s+l 5), was added to G^ to reshape it. 
is shown on a Nichols Chart in Figure 9 and t.._ requency 
plot of the close-loop transfer function, T * , where 


*G # 


). is 


T ;.=^/(l + L;,) = (^ r G;j/(G 

shown in Figure 10 which shows the close-loop frequency 
response of T* after reshaping but with no prefilter applied 

yet. Sufficient gain is available here. A lead compensator of 
(S+l)/(S+2) is added to haunch up the severely deteriorated 
curve at frequencies over 1 rad /sec, and a lag compensator 
of (S-k 25)/(S+.15) is added to steeper the gain curve at low 
frequencies to provide a smoother k/ s curve for good pilot 
handling qualities. The close-loop response after adding the 
prefilter is similar to Figure 8, The prefilter selected is 
15(S+25)(S+l)/((S+.15)(S+2)). 


Results and Discussion 


For bank angle control using only the throttles, p 
feedback was found effective in increasing Dutch-roll 
damping. Bank angle feedback is crudal to lateral phugoid 
damping. Yaw rate feedback, on the other hand, increases 
Dutch-roll damping very little and actually decreases the 
lateral phugoid damping. Tables 4 and 5 compare the 
dynamic modes of the bare airframe with those from the 
augmented control implemented in a previous simulation 
study (References 1 and 2 ), and with those from the QFT 
implementation. 


Table 3. Longitudinal Mode Comparison 



Density 

Phugoid 

Short 

Period 

Engine 

Pre- 

Filter 

G e “' 


K y 

K , 

Bare 

Airframe 

(1.4E-6) 

(.04..13) 

(.65,1.4) 

(.55X5.2) 

* 

m 

X 

X 

m 

Simulation 

Augmented 

Control 

(4.7E-6) 

(.52,-24) 

(52,13) 

(.4X5.2) 

10 

1 

10 

1 

4 

QFT 

Augmented 

Control 

(3.4E-6) 

(.62,32) 

(.46,1.6) 

(.3X53) 

6 

16 

1 

1 

1 

j 

60 
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Table 4 Lateral Mode Comparison 



Spiral 

Dutch 

Roll 

Roll 

Engine 

Pre- 

Filter 

<4 

' G ‘ 

K c 

K P 

k p 

Bare 

Airframe 

a.iE-4) 

(.12,-99) 

(1) 

(.55X5) 

* 

* 

* 

>fr 

* 

* 

Simulation 

Augmented 

Control 

(.73,35) 

U5,.99) 

(1) 

(5) 

40 

1 

1 

.5 

1 

3 

QFT 

Augmented 

Control 

(.39) 

(.29,1.0) 

(13) 

(.45X5) 

2.5(.25) 

(125) 

1 

1 

LI52 

(1.5) 

4 

* 


For longitudinal control, pure gain compensation 
was used. The short period mode has a frequency near 15 
rad /sec, which is beyond the frequency that the throttles 
can control. Therefore, the primary concern was to increase 
phugoid damping and frequency. The short period 
damping decreased from .52 to .46 rad/sec while the 
phugoid damping and frequency increased from .52 to .62 
and from .24 rad/sec to .32 rad /sec, respectively. This 
increase of response frequency can also be easily seen on the 
flight path angle response (not shown here). 

For laTSral control where pole/zero compensation 
was used, the Dutch-roll damping was almost doubled, 
from .15 to .29. The "simulation augmented controller" 
caused a lateral phugoid mode, (.73, .35], which combines 
the spiral and the slow engine mode. The QFT controller 
eliminated the lateral phugoid mode and resulted in higher 
damping for the Dutch roll mode (0.29 versus 0.15). 

Turbulence Response 

The response of the QFT flight path controller 
under intermediate turbulence was excellent. Since gusts 
were input, more than one simulation run was made to 
examine the tracking integrity under turbulence. Both 
controllers performed well, but the lateral QFT controller 
showed undue sensitivity to Kp during bank angle tracking. 

During investigation of the bank angle tracking 
problem, it was found that the p being feedback in the 
simulation was the p at the c.g. instead of at the nose 
boom. The nose boom P was then modeled into the B-720 
simulator, which improved the bank angle tracking under 
turbulence. The P at the nose boom has two extra terms, one 
a function of roll rate, the other a function of yaw rate. It 
was thought that the extra yaw rate term might have 


stabilized the bank angle tracking (the roll rate term was too 
small and was neglected.) 

System response to configuration variations, for 
flight path control and for bank angle control, was very 
good. The robustness of the flight path control was 
improved by QFT. The Dutch-roll oscillation in the original 
simulation compensation was taken out by QFT 
compensation; however, the tracking did not improve. 

Conclusions and Recommendations 

For throttles-only pitch control using a QFT 
controller, the control bandwidth , tracking and control 
robustness were improved by QFT. For bank angle 
control, QFT has improved the Dutch-roll oscillation 
problem and performed well under no turbulence. 

However, the lateral phugoid tracking under intermediate 
turbulence did not perform well. Apparently a compromise 
is required between Dutch-roll and lateral phugoid mode 
damping. Further investigation is recommended for bank 
angle tracking under turbulence. The impact of system non- 
linearities, such as rate and thrust limits, was significant and 
resulted in a decrease in the bandwidth specification used in 
the QFT analysis. 
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c ?(de*/*x) _ >cc) 

*( lbs) *( ibs) iA kmi 

i (lbs) 2(lbs) 


2(lbs) 


*(lbs) 


236E-04 

(0) 

(-1.17E - 05) 

vo. *0) 

(0.61) 

x) 233E-04 
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(0.635, 

0.563) 

1.976E- 04 

<0) 
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1.955E - 04 
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(0.819 , 

0.508) 


Z796E-Q5 
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(0370 , 

3.008) 

-1.819E- 05 
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(2.255) 
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2.130E - 05 (0.167) 

(0351 , 

3.038) 


1.470E * 05 
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(0.261) (0.460 , 3.426) 


(l 438E - 05) 

(3.918E-02 , 
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1301) 


(1.878E-05) 

(7.190E - 02 , 
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(0588 , 

1279) 


Lateral Transfer Functions: 


i(lbs) ~ A 2<lbs) 
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-158E-03 (-.0805) 
N fi(4cg) m ’ 7-M22) 

x(lbs) * -IVif -*03 (V0723)* 

-liSif L’03*“ 7-6*19)' 


(.92 7) 
7904) 

Tiii) 

7 - 940 ) 


3.l9f-04 ( 468 , 3.65) 

^#(dcg) _ ITi? '- 04 "*<iVr"4.T7) 

i(lbs) *Yg9£-“04“‘ (.44773 96) 

Yo4E~-Q*~~ ( 593 ’ _ 4 33 ) 


(.0001) 

0.01) 

(.116 . 1.05) 

(.006) 

(1.05) 

( 067 , 0.93) 

(.0028) 

(1.06) 

(.114, 1.08) 

(.0065) 

(1.09) 

(.060, .944) 
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